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Abstract 

We demonstrate how to make rigorous predictions for electroweak matrix elements in nuclear 
systems directly from QCD. More precisely, we show how to determine the short-distance con- 
tributions to low-momentum transfer electroweak matrix elements in the two-nucleon sector from 
lattice QCD. In potential model descriptions of multi-nucleon systems, this is equivalent to uniquely 
determining the meson-exchange currents, while in the context of nuclear effective field theory, this 
translates into determining the coefficients of local, gauge-invariant, multi-nucleon-electroweak cur- 
rent operators. The energies of the lowest-lying states of two nucleons on a finite volume lattice 
with periodic boundary conditions in the presence of a background magnetic field are sufficient to 
determine the local four-nucleon operators that contribute to the deuteron magnetic moment and 
to the threshold cross-section of np — > dj. Similarly, the energy- levels of two nucleons immersed in 
a background isovector axial weak field can be used to determine the coefficient of the leading local 
four-nucleon operator contributing to the neutral- and charged-current break-up of the deuteron. 
This is required for the extraction of solar neutrino fluxes at SNO and future neutrino experiments. 
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I. INTRODUCTION 



A major goal of nuclear physics is to be able to rigorously compute the properties and 
interactions of nuclei directly from QCD. While QCD is formulated in terms of quarks 
and gluons, the relevant degrees of freedom of strongly interacting systems probed at low- 
momentum are the lowest-lying hadrons. While the properties of these hadrons are well 
known experimentally, and some of their interactions are also well studied (such as the 
nucleon-nucleon (NN) cross-section), in general, direct measurement of the interactions that 
contribute to nuclear processes is not possible. For instance, there are multi-nucleon interac- 
tions induced at the chiral symmetry breaking scale, that will not be directly measured, but 
which contribute to the interactions of nuclei. The only way of rigorously computing strong 
interaction observables and electroweak observables involving hadrons, is with lattice QCD. 
In lattice QCD space-time is discretized and Green functions are evaluated in Euclidean 
space. At this point in time the variety of processes that can be addressed with lattice QCD 
is quite limited. Present day computational power restricts the sizes of lattices, the lattice 
spacings, and the quark masses that can be used in simulations, making the extraction of 
physical observables non-trivial. 

Even if infinite computing power were to become available tomorrow most of the formal 
framework with which to calculate most nuclear properties and observables does not yet 
exist. Notable exceptions are the NN and nucleon-hyperon (NY) scattering lengths and 
range parameters [J 0, IE 0, HI- At the heart of the issue is the Maiani- Testa theorem [f| 
which precludes determination of scattering amplitudes from Euclidean-space Green func- 
tions at infinite volume away from kinematic thresholds . By generalizinga result from 
non-relativistic quantum mechanics £ij to quantum field theory, Luscher 0,1^ showed how 
to extract the 2 — ► 2 scattering amplitude from the energy-levels of two particles in a finite- 
volume lattice simulation. This technique has been used to determine the low-energy tttt 
phase shifts directly from QCD, e.g. Ref. |7[. However, only one lattice QCD calculation of 
the NN scattering lengths |8[ exists, and it is quenched with relatively large quark masses 
(for a review see Ref. |9(). 

In this work we establish a framework with which to extract the electroweak matrix ele- 
ments in the two-nucleon sector. This is only a small step toward a framework with which 
to determine electroweak matrix elements in arbitrary multi-nucleon systems, however, the 
two-nucleon sector is pleasantly challenging all by itself. Naively this appears to require 
the computation of hadronic five-point correlators, a daunting task that is beyond the reach 
of current computers. However, calculation of the energy-levels of two nucleons in a finite 
volume in the presence of appropriate background electroweak fields allows for these matrix 
elements to be extracted with little more computing resources than are needed to calculate 
the scattering parameters. We start by studying the energy-levels of two nucleons immersed 
in a background magnetic field. Clearly a calculation of the deuteron electromagnetic prop- 
erties will not be heralded as a major accomplishment (though it would be an imp ressive test 
of lattice QCD) as they are well-known from precise experiments, e.g. Ref. [l0|. However, 
they do provide an illustration of how the long-distance and short- distance contributions 
can be isolated with lattice QCD. This technique is then generalized to weak background 
fields in order to allow for calculation of the axial matrix element that dominates vd break- 
up III) 2 3 , the process used by SNO [l4| to determine the flavor composition of the solar 
neutrino flux. The present-day uncertainty in this matrix element will provide a significant 
uncertainty in the determination of neutrino mass differences and mixing matrix derived 
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from future experiments. 

Determinations of the electroweak matrix elements will be independent of the hadronic 
theory used to compute the infrared (IR) properties of QCD as long as it is a complete 
theory consistent with all the symmetries of QCD. In this work we use effective field theory 
(EFT), and in particular, the pionless EFT, EFT(ft) to study the very low-momentum 
(|p| <C m^) interaction of nucleons and electroweak gauge fields 2, 13, Hf. EFT(yf) has 
been developed over the last few years to describe systems with unnaturally large scattering 
lengths, as are found in the multi-nucleon sector of QCD. Such large scattering lengths result 
from a fine-tuning that nature has presented to us, more than likely for anthropic reasons, 
and as a result, the analysis in this work will not apply for arbitrary values of the light-quark 
masses. 



II. ELECTROMAGNETIC OBSERVABLES IN THE TWO-NUCLEON SECTOR 
AND "MESON-EXCHANGE CURRENTS" 

In the traditional approach to calculating electroweak observables in nuclear physics, one 
assumes that all processes are dominated by one-nucleon interactions with the electroweak 
probe, and that multi-nucleon contributions are small. These multi-nucleon interactions 
result from currents coupling to the exchange of single mesons, 7r's, p's and so forth, the 
same mesons that are used to construct nucleon-nucleon potentials 1 and are called "meson- 
exchange currents" (MEC's) 2 . In the context of EFT, these matrix elements receive contri- 
butions from one-nucleon interactions with the electroweak current, from interactions with 
the mesons that have masses below the cut-off of the theory (in the pionless theory there are 
no mesonic contributions, while in the pionful theory there are contributions from pions) and 
from local, two- (or more) nucleon - electroweak current interactions. These multi-nucleon 
interactions result from the physics at momentum scales above the cut-off of the EFT that 
has been integrated out. A significant advantage of the EFT is that no attempt is made to 
force these interactions to look like they result from the exchange of mesons - they are what 
they are. However, before the EFT can become predictive, the coefficients of the various 
operators that occur in it must be determined either from experiment or from lattice QCD. 

As the scattering lengths in both the and 3Si channels [a\ = —23.714 fm and = 
5.425 fm) are unnaturally large (both channels are near their IR fixed-points |15l ll6lll9ll20l|). 
EFT(^") can even be used to describe electroweak processes involving the deuteron. EFT(^) 
can be recast into a somewhat simpler theory with which to perform computations up to a 
given order in the momentum expansion by introducing dynamical di-baryon fields |2l|, |22j . 
As S-matrix elements calculated in the usual formulation of EFT (ft) and with the di-baryon 
formalism are the same at any given order in the expansion and the two frameworks differ 
only in the ultra-violet (UV), both forms will give identical results for observables in a 
finite- volume. 

In this work we will not reproduce the methodology of EFT(^) |l8j . however, we shall 



1 The modern potentials that best fit all available data, e.g. AVis, do not rely on single meson exchange to 
describe the short-distance component of the NN potential, instead use some "best-fit" functional forms 
that minimize the \ 2 /dof of the fit to NN data. 

2 Of course, this terminology assumes that all interactions can be represented by meson exchange, which 
while true in the large- Nc limit of QCD, is not true in general. 
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discuss relevant aspects of the di-baryon formalism (2ll . 2^ |. In terms of nucleon and di- 
baryon degrees of freedom, the leading order (LO) low-energy strong interactions for |p| <C 
m^/2 are described by a Lagrange density of the form 
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where N is the nucleon annihilation operator, P is the ^Si di-baryon annihilation operator 
with spin-index j, and s a is the ^0 di-baryon annihilation operator with isospin-index a. 
Further, P| is the spin-isospin projector for two nucleons in the channel with spin-index 
j, while Pi is the spin-isospin projector for two nucleons in the So channel with isospin-index 



r 2 ®a 2 a j , Pi 
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(2) 



In the Lagrange density in eq. (0), a factor of the nucleon mass M has been absorbed into 
the definition of the nucleon fields and similarly a factor of 2M into the di-baryon fields. 
The S-wave interactions are enhanced by a factor of the expansion parameter, 1/Q, and 
are treated non-perturbatively. However the interactions that introduce mixing with higher 
partial waves, e.g. Si — 3 Di mixing, are suppressed by at least Q 2 and so only S-wave to 
S-wave interactions are required to the order we are working 3 . In order to recover the 
scattering amplitudes in both S-wave channels, the constants that appear in eq. are 
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where /1 is the renormalization scale, 03 and r 3 are the scattering length and effective range 
in the channel, and a\ and r\ are the scattering length and effective range in the X S§ 
channel. 

For the processes under consideration, the local, gauge-invariant, electroweak interactions 
between two-nucleons and the gauge fields need to be included in addition to the interactions 
that result from gauging the Lagrange density in eq. The additional electromagnetic 
gauge-invariant operators that contribute to the processes of interest at LO and NLO in the 
power-counting are jzj |iF 
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where k p = k,q + K\ and k„ = k,q — k± are the proton and neutron magnetic moments, e is 
the electric charge, and B is the magnetic field. The magnetic moment interaction of the 



The analysis would become far more involved if this suppression was not present - see Ref. [37 
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di-baryon is written in such a way that if Z 2 = the deuteron magnetic moment is the 
sum of the neutron and proton magnetic moments (HI, 21. 

HUH (at NLO) 



^ d = T7 ( K ° + — ^ ) ' ( 5 ) 
M \ 1 - 7 r 3 J 

where 70 is the binding momentum of the deuteron in an infinite volume with no background 
fields, and satisfies the truncated effective range expansion, 

1 1 2 

7o 7/37o = , (6) 

where the deuteron energy is E = — Jq/M. As spin-dependent interactions are suppressed 
in the infrared, I2 is expected to be small and this is verified experimentally. To make 
some connection with traditional nuclear physics, the constants I2 and l\ correspond to 
multi-nucleon interactions with the electromagnetic field induced by MEC's. 

The breakdown scale of EFT(tF) restricts the magnitude of electroweak fields that can be 
described. In order to have a convergent expansion, the magnetic field must be less than 
|eB| < 2m n M ~ 10 5 MeV 2 so that the one-nucleon magnetic moment interaction shifts the 
nucleon mass by less than m n . Further, requiring that the nucleon magnetic polarizability 
shifts the nucleon mass by less than the magnetic moment (see eq. (J7J) below) requires that 
Aix^Bl < fi p , n B . Using the experimental values of (3 P = 1.2 ± 0.7 ± 0.3 ± 0.4 x 10" 4 fm 3 
and /3 n = 6.5±2.4±3.0 x 10" 4 fm 3 0,113 give a constraint of |eB| < 1.8 x 10 5 MeV 2 , which 
is a somewhat weaker constraint than from the magnetic moment interaction. A more naive, 
but still reasonable restriction is that |eB| <C m 2 ~ 2 x 10 4 MeV 2 arrived at by requiring 
nothing more than that "|eB| is small", and it is this more restrictive limit that we adhere 
to. 



III. LATTICE QCD AND BACKGROUND FIELDS 

A direct calculation of multi-nucleon processes such as np — > dj is impossible on a 
Euclidean space lattice. As briefly mentioned in the introduction, the obstruction is the 
Maiani- Testa theorem which precludes the calculation of S-matrix elements away from 
kinematic thresholds in an infinite volume. In principle however, by performing lattice cal- 
culations of matrix elements of the relevant electroweak operators at unphysical kinematics 
in appropriate multi-nucleon states, e.g., (d(p = 0)| J^\n(p = 0)p(p = 0)), the undetermined 
coefficients of the low-energy EFT (modified to accommodate the injection of energy by the 
lattice) can be determined, leading to an ab initio calculation of the low-energy dependence 
of the process in question. Unfortunately, by present day standards, calculation of such 
five-point correlators is prohibitively expensive in computational terms. In order to do such 
a calculation, the two-nucleon (six quark) sources and sinks and the electroweak current 
operator must be tied together with quark propagators in all possible ways and the corre- 
sponding contributions averaged over the ensemble of gauge configurations. Naively, this is 
one or two orders of magnitude more difficult than computing the NN four-point correlator. 
Clearly, a more cost-effective solution to this problem is desired 

We will show that calculations of electroweak processes, such as np — > and v e d — > 
nne + , will be feasible in the near future. Our method is to some extent motivated by the 
earliest lattice calculations of the magnetic moments of the proton and neutron. Modern 
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calculations of these quantities are performed by taking nucleon matrix elements of the 
electromagnetic current operator. However, in the early days of lattice QCD, calculations 
of nucleon matrix elements were computationally prohibitive (just as multi-nucleon matrix 
elements are today) and the magnetic moments were first extracted by measuring the shift 
in the nucleon mass when the lattice is immersed in a background magnetic field 2y, 27. 28|. 



Taking the magnetic field to be constant and aligned along the z-direction, B = B e z the 
nucleon masses are shifted, e.g. the mass of a proton with spin aligned in the direction of 
the field is 

M^(B ) = M + ^ p B + A7r^B 2 + O(B 3 ) . (7) 

However, it is the energy of the nucleon that is measured on the lattice, and when immersed 
in a background magnetic field the energy eigenstates of the proton in an infinite volume 
are 

E$\Bo) = M + + ^ + » P B + 4vr|i? 2 + O(5 3 ) , (8) 

where n is an integer denoting the Landau-level occupied by the proton, and p z is the 
momentum of the proton parallel to the magnetic field. For weak magnetic fields, the spin- 
dependent energy-splitting is proportional to the magnetic moment of the proton, u p , while 
for stronger fields the magnetic polarizability, (3 P becomes significant and measurable |29, 3Q| . 
To be complete we note that by using background electric fields, electric polarizabilities 
31 . 32l ] and the electric dipole moment of the nucleon 33j have been investigated. 



To apply a spatially constant background magnetic field aligned in the z-direction in a 
lattice simulation, one takes a SU(3) gauge-configuration and modifies the link fields by 

U^x) -> U^x)U^(x) , (9) 

with 

U$(x) = 1 U{ xt (x) = e +lf3x2 U* xt {x) = e~ i(3xi , (10) 

where (3 is proportional to eB Q b 2 where b is the lattice spacing 4 . A vector potential with 
periodic boundary conditions in the transverse directions requires that eB A xy /27T G Z, 
where A xy is the transverse area of the box. This requirement generates large magnetic 



fields for present day lattices j2y, |27|, |28| , large enough to shift the nucleon mass by > m^, 
and thus too strong to be described by EFT (7^. Not requiring periodic boundary conditions 
on the vector potential introduces an exceptional pla que tte that can be placed wherever it 
does the least "harm" in the single nucleon case |2y, |27|, |28(. It remains to be seen if this 
plaquette can be placed somewhere to have minimal impact on multi-nucleon simulations or 
if the periodic requirement must be enforced. 

For the weak interactions, and in particular for an axial-vector interaction, the left- 
and right-handed quarks need to be immersed in different background fields. One way to 
implement this is to use domain- wall (Kaplan j34|) fermions in the simulation and have the 
background field vary only in the fifth dimension. A background axial field was used in the 
first lattice calculation of the axial coupling of the nucleon |35| . 

By applying such background electroweak fields to nucleons in a finite volume, the energy- 
levels become sensitive to the various one- and two-nucleon electroweak current operators 



We work in the continuum limit, b — > 0, throughout. 
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that contribute to the deuteron properties and interactions, such as np — > d'y, v e d — > nne + 
and other related electroweak processes. It is interesting to note that lattice calculations of 
the energy-levels relevant to the deuteron magnetic (or weak) moment (isoscalar operator) 
will be more difficult than those impacting isovector processes as they require the evaluation 
of disconnected diagrams (where the operator does not connect to valence quarks) which 
are significantly harder to compute, even in the single nucleon sector |36| . 

As we will discuss later, it will be necessary to work with asymmetric volumes of dimen- 
sions r\\L ® 772L (g) L (3?}. The energy-levels of two-particles in an asymmetric volume have 
been discussed recently j3?J and we review this analysis and extend a number of relevant 
aspects in Appendix [XJ It is straightforward to show that the locations of the low-lying 
energy-levels of two-nucleons interacting in an S-wave on an asymmetric lattice are given by 
the solutions of H 0, H 

pcotS- ^5(r/!,r/ 2 ;p 2 ) = , (11) 

where p cot 5 uniquely describes the infinite- volume S-wave scattering amplitude below in- 
elastic thresholds. The function S(r]i,r]2]p 2 ) is given by 



1 An 1 

S(vi,V2;p 2 ) = ^ - 4vrA n , (12) 

V1V2 — l n l — P 



where | n| 2 = ^-n^ + ^-n 2 , + n 2 , and the limits of the linearly-divergent, three-dimensional 
sum are given by the ellipsoid |n| 2 < A 2 . The energy corresponding to a given value of p is 

E= M =V MV ■ (13) 

Given the behavior of two-nucleon energy-levels as a function of the volume, the phase shifts 
(or at least the first few parameters of the effective range expansion p cot 5 = — - + + . . .) 
can be recovered. To aid in such recovery, one can find analytic, large volume (L — > 00) 
expansions for the locations of the energy-eigenstates, as done by Liischer for the case of 
771,2 = 1 [HIH El IE H[ and done by Li and Liu j33] for the lowest continuum level in asymmetric 
volumes. We detail this construction in Appendix |X[ and give explicit expressions for the 
large volume locations of the low-lying levels in the two-nucleon sector. A small L expansion 
of eq. (fTTj) can also be performed |5| , though we do not do this here. 

It is straightforward to extend this analysis of two-nucleon energy-levels to include a 
background electroweak field 5 . The presence of the background fields can give rise to a 
number of complications such as Landau-levels and mixing of two-nucleon states of differing 
total spin and isospin. For this reason we shall develop the necessary generalizations of the 
above formulas as we need them. 



5 Recently, Bedaque [iisl has considered simulating the NN system in a background field that couples to 
baryon number, but which has vanishing field-strength. It is found that this is equivalent to performing 
the simulations with twisted boundary conditions on the quark fields, and that the energy-levels of the 
system are sensitive to the phase difference between the boundaries. This appears to be a promising 
avenue for future investigations. 
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The Hamiltonian of a structure-less point particle of charge e, moving in a uniform, 
time-independent magnetic field in the ^-direction is 



where to = | | , and l z is the angular momentum operator parallel to the magnetic field. In 
an infinite volume, the energy-eigenvalues are given by the expression in eq. (jHJ)- However, 
in a finite volume the oscillator potential is bounded above, and for small enough lattices the 
potential arising from the magnetic field can be treated as a perturbation. Consequently, for 
a small magnetic field we can work with momentum eigenstates in the transverse direction. 
Requiring that the energy-shifts due to these magnetic interactions are small compared with 
the inter-level spacing of the two-nucleon momentum eigenstates in the finite volume gives, 

\eB \ < — — , (15) 
L x 

where Lj_ is the size of the transverse dimensions. In what follows we work with momentum 
eigenstates in the transverse directions, and thus are constrained by eq. . 



IV. DEUTERON MAGNETIC MOMENT : ISOLATING THE TWO-NUCLEON 
CONTRIBUTION 

It is useful to begin with the simplest case to analyze and understand — the deuteron 
magnetic moment. To isolate the two-nucleon contribution from the one-nucleon contribu- 
tion (the magnetic moments of the proton and neutron combined in the S = 1 channel), 
the energy-levels of two nucleons interacting in the ^S*i channel in a finite volume with a 
background magnetic field in the ^-direction, B = B e 2 , are computed. 

It is obvious from eq. (J7|) that the m = ±1 magnetic sub-states of the NN system will 
have a linear dependence on the magnetic field for small fields. Indeed, the effect of the 
background magnetic field on the single-nucleon states is to shift their masses from M to 

gB gB 
MpT = M ~2I^ Kp ' Mpi = M+ 2I^ Kp 

MnT = M ~^ Kn ' Mnl = M+ ^ Kn ' (16) 

where the proton and neutron magnetic moments. These mass-shifts appear as 

residual mass-terms in the single nucleon propagators that contribute to the infinite- volume 
NN scattering amplitude, and hence to the position of the energy-levels of two-nucleons in 
a finite volume 6 . In the weak-i? limit contributions to the NN energy-eigenvalues that are 
quadratic and higher-order in B Q , such as the nucleon polarizabilities and modifications aris- 
ing from the eigenstates in the transverse directions being Landau-levels and not momentum 
plane-waves, can be neglected. 



6 Only the free-space nucleon mass M is removed by a phase redefinition of the nucleon field in defining the 
non-relativistic theory. The presence of interactions between different spin and isospin components of the 
nucleon field prevents a rescaling of each component to completely remove the mass-shifts in eq. 
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Since it is a non-relativistic system, a simple recipe can be followed to determine the 
location of the low-lying energy-eigenstates of two nucleons in a background magnetic field 
in a finite volume 0, 11 • Construct the inverse scattering amplitude in the background mag- 
netic field, replace the continuum momentum-space loop integral (arising from the bubble- 
diagrams) by a discrete sum over the allowed momenta in the volume and then find the zeros 
of the real part 0, 0. It follows directly from this recipe that the location of the m = ±1 
energy-levels are given by the solutions to 



pcotS 3 - — 5(771, 772; p 2 ±ujj) =F 

TlL 







(17) 



where « = |(k p + K n ) and pcot5 3 (p 2 ) = — — + ^p 2 + . . . is the effective range expansion 
in the 3 Si channel. The one-nucleon contribution to the location of the NN energy-levels is 
through the term 



~2 



u, 



An 2 



eB K 



(18) 



and through k in the last term in eq. (|17|) . The one-nucleon and two-nucleon (determined by 
I2) terms enter eq. (fTTj) differently. Consequently, the short-distance and the long-distance 
interactions with a magnetic field can be separated, i.e. I2 can be isolated from k,q on the 
lattice. Note that the expression in eq. (fT7|) reduces to the expression in eq. (fTT]) in the limit 
of vanishing magnetic field. 

For symmetric lattices with 77^2 = 1, it follows directly from eq. (|17|) that for L ^> a 3 
there are two states that become the deuteron m = ±1 states in the infinite-volume limit, 
that have energies 
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(19) 



in a finite- volume. The deuteron magnetic moment at NLO in the EFT is given in eq. 
but eq. (|19|) is true at all orders in the expansion. The leading finite-volume corrections in 
eq. (fT^j) are the same as those at B = |5[ . Although this gives a clean lattice measurement 
of the deuteron magnetic moment (itself a very good test of the lattice approach), it does 
not facilitate a separation between the one- and two-nucleon contributions. 

Making an expansion of eq. ()17|) about the infinite- volume limit, E = =p eB ^° , the energies 
of the first pair of continuum states in the limit where L>a 3 are 
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where an effective scattering length is defined to be 

1 



a 3 2 



(20) 



(21) 



and the constants c 1>2 for 77^2 = 1 are given in Table ITD 0, 0|. The energies of the second 
pair of continuum states in the limit where L>a 3 are 
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where 



? 1 1 4tt 2 

pcot5 = ~A~ 3 + 2 r3 l? ' ^ ' 

the constants c[ 2 for r] 12 = 1 are given in Table HP 0,0], and A 3 is defined in eq. (J2T]) . 

The energies in eqs. ffeUj) and (j22J) closely resemble the well-known Liischer formulas 0, HI , 
with some interesting modifications. The levels receive shifts proportional to k correspond- 
ing to the shifts that two non-interacting nucleons would receive in a magnetic field in an 
infinite volume. In addition, instead of the physical scattering length entering into the 1/L- 
expansion, it is the effective scattering length, A% that appears. The effective scattering 
length, defined in eq. (}2"T]) . depends upon the actual scattering length and the short-distance 
contribution to the NN interaction due to the background magnetic field. These expressions 
are not entirely unexpected. In the continuum states the nucleons do not "see each other" 
much, and the short-distance interactions are a perturbation (as is clear from eqs. (J2"U|) and 
(J220- The single particle energies dominate and the short-distance interaction, which is 
the sum of both the purely strong interaction and that induced by the background field, 
is inserted perturbatively. Measuring the energies of the continuum levels over a range of 
lattice volumes will allow for a clean separation between the one-nucleon and two-nucleon 
contributions to the deuteron magnetic moment. 

Given that the ^S"i channel is very close to its non-trivial fixed point (03 — > 00), most 
choices of magnetic field will move the system away from the fixed point, and hence \A%\ 
will typically (i.e., for natural values of I2, \fa\ ~ 1 fai) be dramatically smaller than \as\. 



This suggests that, in principle, lattices even smaller than those discussed in Ref. |5| could 
be used to study the two-nucleon sector when a background magnetic field is present. How- 
ever, nature is not as kind as it might have been, suppressing the isoscalar spin-dependent 
interactions (proportional to I2) at low-momentum by even more than the naive l/iV ( H~ J ' 
obtained in the large- Nc limit j39|. The experimental value of the deuteron magnetic mo- 
ment Kd = 0.85741 is very close to the single particle sum of 2k = 0.87976 which results 
in an unnaturally small value of I2 = —0.0576 fm in eq. (JHJ). Therefore, for magnetic fields 
that allow for an analysis of the two-nucleon sector within EFT(^), the effective scattering 
length, A3, is only moderately different from the actual scattering length, as. Specifically, 
for eB = 10 4 MeV 2 , A3 differs from a 3 by 10% (for the physical value of I2), an amount 
that should produce observable shifts in the continuum two-particle energy levels calculated 
in lattice simulations. In fig. ^the energy-levels of the three lowest-lying m = ±1 states 
are shown, along with the analytic expressions in eqs. (J19)) . (J20|) . (J22j) and their asymmetric 
generalizations that are valid for L ^> a 3 . 

An interesting "twist" could be applied to the systems we are considering. It is possible 
that the light-quark masses, m q , could be tuned to values for which the physical scattering 
lengths in the ^S*i or &o channels, or both, become infinite ji^, 41, 4^, 43|. As both channels 



have unnaturally large scattering lengths at the physical values of m q , only small changes in 
these masses will be required, and hence the values of the operators in the pionless theory 
describing the unphysical system, will differ from their physical values by only small amount. 
An advantage of considering such an unphysical system is that A 3 will depend only on l 2 and 
its determination on the lattice would provide a measurement of I2 for the physical theory 
with a reasonably small systematic uncertainty. 
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FIG. 1: The energies of the lowest-lying %i states with m = ±1 on the lattice as a function of the 
lattice size, L, in the presence of a background magnetic field of strength cBq = 10 2 MeV 2 . The 
left plot is for a symmetric volume, while the right plot is for r/i = 7]2 = 0.2. The experimental 
values of the nucleon and deuteron magnetic moments have been used. The two types of triangles 
correspond to numerical solution of eq. \17\ ), and the curves correspond to the large-L asymptotic 
behavior given by eqs. hiy\) . \2U\) and h22\) and their asymmetric generalizations. The gray shaded 
regions denote lattice sizes that are large enough so that the effects of the transverse eigenstates 
being Landau-levels and not momentum eigenstates cannot be ignored. 



V. 



np 



dj : ISOLATING THE TWO-NUCLEON CONTRIBUTION 



The radiative capture process np —>■ d*y, and its inverse reaction (photo-disintegration 
of the deuteron) provide a classic demonstration of the importance of MEC's in nuclear 
physics. The one-nucleon contribution to thermal neutron capture provides only ~ 90% of 
the experimentally measured cross section. From the EFT perspective, it is natural that 
two-nucleon electromagnetic interactions can contribute to this process, in the same way 
that they contribute to the deuteron magnetic moment. The cross sections for np — > d'y 
for cold neutrons and for the photo-disintegration of the deuteron near threshold are both 
sensitive to 1%. The radiative-capture cross section near threshold is given by [22 



a(np — > d^) 



Ana (70 + 


IP 


i 2 ) 3 


M 4 7o 3 | 


Pi 





I^Mll 2 + |^£l| 2 + 



(24) 



where contributions from multipoles other than Ml or El, denoted by the ellipses, are para- 
metrically suppressed at low-energies. Near threshold, this process takes the two-nucleons 
between the &g channel and the channel and the cross section is dominated by the Ml 
amplitude jij" 



X 



1 



1 



Ml 



1 119 

Fi|p| 2 



1 p 



«i7o 



7o + IPl 



1 1 , 

7o + 7T r i P 

ai 2 



+ 7 7 ° 

+ /l T 



.(25) 



It is worth stressing that this very s imp le expression for the cross section (including the 
analytic form for the El amplitude |44j, |45| ) describes all available low-energy data with 
high precision 1(1 HI, 45^ . 
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A background magnetic field aligned in the z-direction induces a mixing between the l So, 
I z — state and %>i, m = state. Therefore, the energy-eigenstates in the background field 
will be a mixture of both the ^So an d the channels, which gives rise to quite interesting 
behavior. In large volumes, the continuum states, for which the strong interactions play only 
a small role, are perturbatively close to the single particle eigenstates anc l n^p^, with 
energies dictated by the single particle interactions with the background field, i.e. linear 
in B Q for small B . Hence, the energies of the continuum states will be insensitive to the 
value of l\. However, the energies of the near threshold bound-states, the deuteron and 
the virtual bound di-nucleon state in the 1 Sq channel, for which strong interactions play an 
important role, will depend quadratically on Bq for small Bq and will have sensitivity to 
l\. This results from the fact that these states have equal admixtures of the single particle 
eigenstates in infinite volumes. Thus, it is clear that we will need to observe energy-shifts 
that are quadratic in B in order to extract l\. This presents a potential problem that we 
did not encounter for the energy shifts of the m = ±1 states. In the directions transverse to 
B, Landau- levels are the single-particle eigenstates rather than momentum plane- waves. As 
discussed in section 11111 this is a complication in infinite volumes, but becomes a far more 
serious issue in finite volumes for systems with H(3) symmetry. In order to avoid having to 
contend with this problem, we eliminate it; highly asymmetric lattices, for which 771,2 ^ 1> 
so as to satisfy eq. (|T5|). strongly suppress excitations in the transverse direction, leaving the 
system dominated by fluctuations in the z-direction. Ideally, one would take the transverse 
dimensions almost to zero to eliminate these contributions almost entirely. However, the 
pion mass sets a lower-limit of rji^L ~ 2 fm in order for thepionless theory to provide a 
valid description. On a symmetric lattice, as studied in Refs. [H 0, 0, 0, S| , the lowest-lying 
states are in the A\ representation of H(3) and as such their energy is dependent upon 
NN interactions with J = 0,4,6,8,.. |46[ with the notable omission of interactions with 
J = 2. On an asymmetric lattice where H(3) is not a symmetry, NN interactions with 
even angular momentum will contribute. However, in the power-counting of EFT(^), higher 
angular momentum contributions are suppressed and can be systematically included when 
required. 

As the background field induces mixing between the and 3 Si channels, the energy- 
eigenvalues are found by diagonalizing the coupled channels system. It is straightforward to 
show that these eigenvalues are solutions to 



p cot Si 



S± + S2 
2vrL 



P COt £3 — 



Si + S2 
2ttL 



eB li Si — S2 



2ttL 



(26) 



where 



Si = S(r]i,r] 2 ;p +u ± ) 



So 



S(vi,V2;P -Ui) 



r.2 



U 



4tt 



eB Ki 



(27) 



K i = ( K p ~ K «)/2 and pcot Si = — + yp 2 + ... is the & effective range expansion. The Ui 
contributions result from the one-nucleon interactions with the background field, while the 
two-nucleon-background field interactions are described by 1%. 

The L 3> 0,1,3, r i,3 asymptotic region is more complicated in these mixed channels than in 
a single channel and we do not have practicable analytic expressions for the energies of the 
lowest-lying levels in this limit. Thus, all of our intuition and understanding has come about 
through numerical solution of eq. ()26|). In fig. El we show the spectra of energy-eigenvalues 
of the %(I Z = 0) - 3 Si(m = 0) NN system in magnetic fields of eB = 4 x 10 3 MeV 2 
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FIG. 2: T/ie energy of the lowest-lying states in the coupled ^3q(I z = 0) — 3 5i(m = 0) channels 
on an asymmetric lattice with r]i2 = 0.1 as a function of the lattice size, L, in the presence of 
a background magnetic field. The left panel shows the spectrum for cBq = 4 x 10 3 MeV 2 , while 
the right panel shows the spectrum for cBq = 5 x 10 2 MeV 2 . The shaded regions indicate energies 
outside of the EFT(yf) and states therein should be treated with extreme caution. 



and eBo = 5 x 10 2 MeV 2 on an asymmetric lattice with 771,2 = 0.1 as a function of lattice 
size, L. As a function of B and h, the continuum levels behave precisely as one would 
expect, and are very insensitive to l±. However, as anticipated, there are one or two states 
in the theory (depending on the choice of parameters) which are very sensitive to l\. The 
location of these states in an infinite volume are finely tuned (near the IR fixed point), and 
as such, "squeezing them" into a volume that is small in the transverse directions destroys 
this fine-tuning. However, the background field, through the operator with coefficient 1%, can 
restore the fine-tuning and make the location of the state depend strongly upon l\. Further, 
the observation of two bound-states in the spectrum that are sensitive to l\ for some values 
of background field and for some volumes indicates that the virtual bound state in the x Sq 
channel can be moved onto the first sheet by the background magnetic field. In the infinite 
volume limit one of the /i-dependent states becomes the physical deuteron. However, in a 
finite volume and in the background field its composition, and the composition of the other 
states, will be quite complicated, each being a combination of the "deuteron", the ^-"di- 
nucleon',' and the continuum states, as deduced from fig. 121 As EFT(^") is applicable only 
at very low-energies, the energy-spectra in fig. |2]are reliable only for \E\< 10 MeV. 

For eBo = 4 x 10 3 MeV 2 , as shown in the left panel of fig. there is one "reliable" 
state 7 that moves to higher energy as a function of increasing volume, and is very sensitive 
to l\ in the physical range of l\. For lattice sizes between ~ 30 fm and ~ 50 fm in the 
longitudinal direction, this level is one of the low-lying ones in the spectrum, and therefore, 
somewhat easier to extract from lattice calculations, than if it were high in the spectrum. 
Also present in this spectrum is the curious feature of "level pair-production" near the cut- 
off of the theory. In the absence of a background field there are two solutions to eq. (jUJ), 
one corresponding to the deuteron and one to an unphysical state lying outside the theory, 



7 By reliable we mean that it lies within the range of validity of EFT(/f). 
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nonetheless both solutions exist, and the one outside the theory is simply discarded as 
unphysical. The same line of reasoning applies to the deep states in fig. We interpret 
these states as the unbound di-nucleon state in the 1 Sq channel mixing with continuum states 
and the deuteron through the background magnetic field, bringing an additional state down 
from the cut-off into the spectrum. Its sensitivity to l\ suggests that it is a strongly correlated 
NN pair, as opposed to predominantly continuum in nature. For eB = 5 x 10 2 MeV 2 , as 
shown in the right panel of fig. 121 the two lowest-lying states in the spectrum are both 
sensitive to l\. It might well be the case that a precise measurement of the lowest-lying state 
in the spectrum at weaker fields is the easiest way to determine l\ from lattice QCD. 

A general statement can be made about the NN system in background magnetic field. 
There are values of the magnetic field, lattice asymmetry and lattice size for which the 
spectrum of energy-eigenvalues allows for the short-distance and long-distance contributions 
to electromagnetic matrix elements to be separated, as is clear from fig. |21 While a lattice 
extraction of these quantities will be very important to perform and provide a proof in 
principle of this method, it will not, unfortunately, tell us anything we didn't already know 
experimentally. 

VI. WEAK INTERACTION OBSERVABLES IN THE TWO-NUCLEON SECTOR 

Of central importance to present day and future neutrino experiments is the ability to 
separate the flavors of neutrinos entering the detectors. The method to accomplish this 
implemented by the SNO collaboration is to use the weak-disintegration of the deuteron, 
e.g. v x d — > v x np [x = e,fi,r) and V e d — > nne + . These processes "look" like the photo- 
disintegration of the deuteron discussed in the previous section when the electromagnetic 
field is replaced by a weak gauge field. There are both one-nucleon and two-nucleon con- 
tributions to these processes, but unlike the electromagnetic case, the two-nucleon contri- 
butions have significant uncertainties associated (for a nice discussion see Ref. [l]]|) with 
them that will impact the future determinations of the neutrino mass and mixing matrices. 
Therefore a determination of the short-distance contributions to the matrix elements for 
weak-disintegration of the deuteron from lattice QCD is important. 

The procedure for determining the two-nucleon weak interactions is similar to that de- 
scribed in the electromagnetic sector with a few minor differences. For the electromagnetic 
interactions it sufficed to use a background magnetic field. However, for the weak interac- 
tions it is convenient to introduce a fictitious SU(2)l® SU(2)r weak gauge symmetry under 
which the up- and down-quarks transform, giving rise to the interactions 

£ int ' = -i [ g L Wl L (uj^l - l5 )u - Sy (1 - 7B )d) 

+ QbWIr («7 m (1 + 7s)« - dY(l + 7s)d) ] , (28) 

and give vevs to W% L and W^ R , such that {qlV^zl) = —(9rW^ r ) = gW. Therefore, in this 
background weak field there is an additional contribution to the strong interaction Lagrange 
density of 

5C = - l -gW (u7* 7fi u-d7*75<0 • (29) 
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This interaction leads to interactions in the di-baryon formulation of FFT(fi) of the form 



5C 



gW — iVVr 3 iV - 



gWl 



l. A 



tls 3 + h.c. 



+ 



(30) 



where the ellipses denotes terms higher order in EFT(yf), and = 1-26 is the nucleon axial 
coupling constant. We note that due to the nature of the interaction, weak "Landau-levels" 
are not present in this system. 

The cross sections for the weak-disintegration of the deuteron and other two-nucleon 
weak processes depend upon the dimensionless coefficient 1\,a- These processes have been 
studied extensively by Butler and Chen jij IHf and also by Butler, Chen and Kong ji^ and 
our (renormalization-group invariant) 1\a is related to the constants they introduce. The 
coefficient of the four-nucleon operator in EFT(/f), ^L^a, is related to I^a via 



2(^-7) 
Cf 0) (/i) L 



r 3 



(31) 



where C^°V) 

are the coefficients of the zero- and two-derivative strong interaction opera- 
tors in EFT(^). Numerically, eq. (J3T]) reduces to 



l lA = -13.4 + 0.27 *Li 



(32) 



when the RG-scale \x = is chosen, and where ^L\ )A is in units of fm 3 . Calculations 
have been done in which a chiral expansion of the weak currents are performed in a manner 
consistent with Weinberg's power-counting [Hoi] . Matrix elements of these operators are taken 
between wave-functions generated with the best modern NN potentials [l3j]. This method is 
somewhat ad hoc, but it has been demonstrated to be convergent where it has been tested. A 
calculation of tritium /3-decay in such a framework leads to [l3j] *Li i A(m n ) = +4.2 + 0.1 fm 3 . 
Chen, Heeger and Robertson ll( have surveyed the experimental and theoretical constraints 
on ^L l j 4(m 7r ) and have also provided a model- independent determination of this quantity 
from the SNO and Super-Kamiokande data, finding ^L\ yA = +4.0 ± 6.3 fm 3 . 

The cross sections for the processes vd — > unp, Vd —>■ Vnp, and Vd — > e + nn are all 



given in terms of ^-Zvi,a("v)) i n Refe. 12 > HM 4^ - and analytic expressions for each can be 
found there. As these expressions are quite complicated, we do not reproduce them here. 
Numerical values of the cross sections for these processes can be found in Refs. 42, 48i Efll]. 
As an example, at E vV 



10 MeV the cross sections are, in units of 10 42 cm 2 , 



a(u x d -> v x np) 
o{y e d — > e~pp) 



1.76 + 0.056 l ljA 
4.07 + 0.12 l 1A 



, o{y x d 
a(V P d - 



-> V x np) = 1.66 + 0.052 Z M 
e + nn) = 1.93 + 0.059 Z M , (33) 

10% level translates into an 



where x = e, /i, r. Thus, a determination of I^a at the 
uncertainty in these cross sections at the few percent level. 

In the presence of a background weak field of the form given in eqs. f!29|) and ([30)1. the 
I z = component of the ]Sq channel mixes with the m = component of the ^Si channel, 
as is the case in the presence of a background magnetic field. The energy-eigenvalues of 
^(iz = 0) — 3 Si(m = 0) NN system in a finite volume with this weak field are solutions to 



pcot^i — 



S1 + S2 
2ttL 



p cot S3 — 



Ox + S2 

2irL 



'gWl 



1,A 



Si — S2 
2irL 



(34) 
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FIG. 3: The energy of the lowest-lying states in the coupled ^3q(I z = 0) — 3 5i(m = 0) channels 
on an asymmetric lattice with rjip = 0.1 as a function of the lattice size, L, in the presence of 
a background weak field. The left panel shows the spectrum in the presence of a background field 
gW = 3 MeV ; while the right panel shows the spectrum for gW = 6 MeV. The shaded regions 
indicate energies outside of the ^FT(f) and states therein should be treated with extreme caution. 



where 

L 2 

Si = S(rn,ri2]p 2 + w\) , S 2 = S(r)i,r)2;p 2 - w\) , w\ = -— gW g A M .(35) 

In analogy with the electromagnetic case, the w\ contributions result from the one-nucleon 
interactions with the background weak field, while the two-nucleon interactions are described 
by As in the electromagnetic case, we have explored the spectrum of this system through 
numerical solution of eq. ([35)1 . The spectra of two-nucleons in background weak fields of 
strengths gW = 3 MeV and gW = 6 MeV are shown in fig. El as a function of L. In 
complete analogy with the electromagnetic case, the continuum levels are insensitive to I^a, 
while there is one level, or two levels for some parameters, that are sensitive to l\ t A due 
to the fine-tuning in the NN sector. For lattices with size L ~ 40 fm and 77^2 = 0.1 the 
state that is sensitive to Ix,a is either the lowest state or one of the lowest states, depending 
upon the choice of field. Over a reasonable range of I^a, the energies of the levels that are 
sensitive to I^a vary by ~ 1 MeV, which should be measurable. 

It is also worth pointing out that an asymmetric lattice is not necessary to study these 
weak interactions. It is the existence of Landau-levels in the case of a magnetic field that 
forced us to consider an asymmetric lattice to suppress excitations in the transverse direc- 
tions. However, for the weak background field in eqs. ()29|) and (|30j) the transverse eigenstates 
are momentum eigenstates, and as such a symmetric lattice is sufficient. 

It is clear from fig. |3] that such lattice simulations of two-nucleons in a background weak 
field will allow for a determination of I^a- We find this to be quite remarkable. As a final 
comment on the weak properties of the two-nucleon sector, one could explore the weak- 
moment of the deuteron. This too will have contributions from the weak-moments of the 
nucleons, and also from short- distance interactions with the neutral-current gauge-field. By 
examining the two-nucleon spectrum in the presence of an isoscalar weak field, the short- 
distance interaction could be isolated from the long-distance contribution and the deuteron 
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weak moment could be determined. This too would be a spectacular accomplishment for 
lattice QCD. 



VII. CONCLUSIONS 

In this work we have shown that lattice simulations of the two-nucleon sector in a fi- 
nite volume and in the presence of background electroweak fields allow for the isolation of 
the short-distance, multi-nucleon electroweak currents, traditionally called "meson-exchange 
currents". In lattice QCD simulations one has many more "knobs" to turn to get at the 
physics one is interested in. In particular, the dimensions of the lattice can be chosen to 
disturb the fine-tuning in the two-nucleon sector in a controlled way to maximize sensitiv- 
ity to the short-distance electroweak interactions. Remarkably, there can be one or two 
levels in the spectrum that are highly sensitive to these short-distance interactions, while 
the continuum states show little sensitivity. Consequently, for some background fields and 
lattice dimensions, the short-distance interactions can be isolated from the long-distance 
interactions. 

One would like to see a higher-order analysis of these finite volume observables in EFT (ji) 
including the effects of the Landau-levels, one- and two-nucleon polarizabilities and other 
relevant contributions. Further, our analysis has been performed exclusively with the pi- 
onless EFT, EFT(/f). While the result are rigorous where applicable, this work needs to 
be extended to a larger range of energies and momenta, which requires solving the pionful 
theory in a finite volume. Progress toward this goal is being made in the zero- and single 
nucleon sectors, but so far little progress has been made in the multi-nucleon sectors. Re- 
cently, Lee, Borasoy and Schaefer |51[ have performed the first lattice simulations of the 
pionful EFT using Weinberg's power-counting, with an eye to a complete numerical solution 
of multi-nucleon systems 8 . Further, encouraging progress has been made toward setting-up 
lattice simulations in EFT(;f) 0, HI . 



The numerical analysis we have presented in this work has been for the physical values 
of the strong interaction parameters, and consequently, for the physical values of the quark 
masses. Lattice simulations of the immediate future will not be performed at the physical 
quark masses, be they quenched, partially-quenched 9 or unquenched. It is possible that the 
scattering lengths in these unphysical simulations will not be unnaturally large, and the 
pionless theory will not be applicable. Therefore, it is more than likely that the pionful 
theory will need to be solved to more precisely understand how the scattering parameters 
depend upon the light quark masses, and to directly extract the interactions that result from 
physics at the cut-off scale of the pionful theory. 

This work provides the first steps toward extracting the electromagnetic and weak prop- 
erties of nuclei from lattice QCD simulations. We encourage lattice practitioners to explore 
this new area and to determine the resources required to perform such calculations. As 
mentioned earlier, even if infinite computing power were presently available, one knows very 



8 The formal problems associated with the chiral expansion in Weinberg's power-counting uncovered in 
Ref. (see also Ref. Q) have not been resolved in Ref. [isij . 

9 Generically, the non- unitary nature of the quenched and partially-quenched versions of QCD is problematic 
when considering two-particle scattering. However, for energies below ghost meson production threshold, 
the partially-quenched theory is unitary and can be analyzed rigorously [5^ . 
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little about what to compute in lattice QCD in order to make predictions about nuclear 
processes. Therefore, we encourage nuclear theorists to extend this work and the work of 
Refs. [U, HI IE IH El so that a rigorous computational framework for nuclear physics is in place 
when the computing power becomes available. 
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APPENDIX A: TWO NUCLEONS ON AN ASYMMETRIC LATTICE 



The energy-levels and formalism associated with two interacting particles on an asym- 
metric lattice of finite volume, have recently been studied |37T |. We use the asymmetric 
lattice to suppress nucleon momentum modes in the directions transverse to the background 
magnetic field in order to extract the electromagnetic counterterm, I2, i.e. to minimize the 
contribution from the N' A ■ ViV and N^A 2 N interactions which generate Landau-levels in 
an infinite volume. 

Explicit evaluation of the function S(r}i, r] 2 ] p 2 ) defined in eqs. (jTTj) and (jT2j) can be ac- 
complished in a variety of ways. Extending the Chowla-Selberg relation [57[ to asymmetric 
lattices, it is easy to show that for x ^ 



—2lT X + 7T 



i m 



3 — 27r|m|a; 



(Al) 



where |m| 2 = ri 2 m 2 + r]2m2 + m 2 . In order to facilitate evaluation of S(r]i,r]2;p 2 ), and also to 
determine the coefficients in the expansion of the energy-levels as a function of lattice size, 
three sums are evaluated for a given lattice asymmetry. The most convergent sum (~ 1/A^ 
in the UV) is 



K, 



m '/-> 



V1V2 ^ |n|' 



(A2) 



A second sum is 
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(A3) 



and it is important to note that J vim is independent of the value of x, and the expression 
in eq. (|A3|) holds for x ^ 0, and also converges as ~ The third sum is 

4vrA r . 
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(A4) 



which is also independent of the value of x, and the expression in eq. (|A4j) holds for x ^ 0. 
The numerical values of these sums are given in Table for different lattice asymmetries 
771 2- I n terms of these p-independent constants, the function S(r]i,r]2]p 2 ) can be written as 



TABLE I: The constants Immi Jmm, K mri2 for different lattice asymmetries 771,772- 



Vl,V2 


I mm 


Jf]m2 


K 


0.1 , 0.1 


206.456 


217.884 


203.475 


0.2 , 0.2 


20.9815 


56.9542 


50.9151 


1 , 1 


-8.91368 


16.5323 


8.40191 
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+ P 2 J, 



l VlV2 



771772 p 2 



V1V2 
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Ifil 



~ 1 
P 



11012 t£o l fi l 4 (|n| 2 -p 2 y 



(A5) 



which is valid for all lattices defined by 7712, and for all p 2 . Accuracy in the 5th or 6th 
significant digit can be obtained by evaluating a few finite sums that require insignificant 
amounts of computer time. An additional check on our numerics is provided by the integral 
representation 



l^l 2 ^ (2p 2 -l) + — EfTT^ 

771772 \Q\ — P 




(A6) 



In the limit that all lattice dimensions are much larger than any strong interaction length 
scale, including the scattering length, approximate formula for the energy-levels can be 
found [J, 0, HI . An expression for the energy of the first continuum state of two particles 
on an asymmetric lattice has been given in this large- volume limit in Ref. 37|, and we add 
to this work by giving expressions for the bound state (the deuteron) and for the second 
continuum level. We find that the bound state energy of two nucleons in the channel is 



7f 
M 



1 + 



7oL(l - 7 r 3 ) 



g -70-L 1 ]_ e -Jom L _|_ \ C -10V2L 
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(A7) 



In the relevant case of 77^2 <C 1, the first term in the round brackets can be neglected and the 
second two dominate the expression (for other cases such as r/1.2 1 the modifications are 
obvious). Physically, the dominance of the transverse dimensions is reasonable for r/1.2 *C 1 
as the contributions to the kinetic energy of the nucleons in the deuteron from these two 
directions are being reduced, moving the deuteron further from its infrared fixed-point. As 



19 



expected for an attractive interaction, the deuteron becomes more deeply bound as the 
lattice volume is reduced. 

The energy of the lowest-lying continuum state, corresponding to the lowest momentum 
mode in all three directions, is 



En 



1 -01(771,7/2) +03(771,7/2) 



+ 



(A8) 



valid in both the ^0 and channels. The values of the coefficients c\ and C2 for various 
values of r/12 are given in Table ITT1 



TABLE II: Values of the constants 01,2(7/1,772) and c\ (m , y?,) 
7/1,7/2- Additional values of 012(771,772) are given in Ref. [37 1. 



for different lattice asymmetries 



771,772 


01(7/1,772) 


02(7/1,7/2) 


4 (7/1,7/2) 


c 2 (7/i,7/ 2 ) 


d 


0.1 , 0.1 


65.7171 


2111.12 


-3.60224 


-52.944 


1 


0.2 , 0.2 


6.67861 


-99.6627 


-2.3242 


0.950737 


1 


0.5 , 0.5 


-3.61168 


6.65945 


0.73400 


0.24642 


1 


1 , 1 


-2.837297 


6.375183 


-0.061367 


-0.354156 


3 



The next continuum level corresponds to non-zero momentum in the longest direction(s), 
e.g. p = (0, 0, ±^r) for 77^2 < 1. For large L, the energy of this level for 77^2 < 1 is 



47r 2 4tan5 
1 ~ ML 2 ~ rj^ML 2 



1 + 



4(771, 772) tan 5 + 4(771, 7/ 2 ) tan 5 + ••] 



(A9) 



where the phase shift is evaluated at |p| 2 = 4n 2 /L 2 , and d is a degeneracy factor, given in 
Table ILTl Again for other choices of 77^2, the corresponding expressions are easily derived. 
Values of the coefficients c[ 2 are given in Table ILTl for some choices of 771,2- The large values 
of some of the Cj's in Table O indicate that it is the "short" direction(s) that are setting the 
scale of the finite- volume corrections, and not the "long" direction(s). 

Finally we note that there are curves in the 7/1-772 plane for which 0^2(771,772) separately 
vanish, as shown in fig. 0] The point on this plot that is most relevant to our analysis 
is 7/1=7/2 = 0.268494 for which ci = 0. Such "magic boxes" may prove to be useful in 
lattice calculations as the asymptotic formulas in eqs. ()A7j) — f[A9|) have better convergence 
properties at smaller L than those associated with arbitrary lattice volumes. Analogous 
contours exist for the 012(771,772) in the 771-772 plane. 
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